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Abstract 

We review the quantum Yang-Mills condensate (YMC) dark energy models. As the effective Yang-Mills Lagrangian is 
completely determined by the quantum field theory, there is no adjustable parameter in the model except the energy scale. In 
this model, the equation-of-state (EOS) of the YMC dark energy, Wy > —1 and Wy < —1, can both be naturally realized. By 
■ studying the evolution of various components in the model, we find that, in the early stage of the universe, dark energy tracked 
' the evolution of the radiation, i.e. Wy 1/3. However, in the late stage, Wy naturally runs to the critical state with Wy = — 1, 
' and the universe transits from matter-dominated into dark energy dominated stage only at recently z ~ 0.3. These characters 
are independent of the choice of the initial condition, and the cosmic coincidence problem is avoided in the models. We also 
I ^ find that, if the possible interaction between YMC and dust matter is considered, the late time attractor solution may exist. 
In this case, the EOS of YMC must evolve from uij, > into Wy < —1, which is slightly suggested by the observations. At the 
same time, the total EOS in the attractor solution is wtot = — 1, the universe being the de Sitter expansion in the late stage, 
and the cosmic big rip is naturally avoided. These features are all independent of the interacting forms. 

PACS numbers: 98.80.-k, 98.80.Es, 04.30.-w, 04.62. -|-v 
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I. INTRODUCTION 



The observations of Type la Supernova (SNIa) [ij, together with the cosmic microwave background radiation 
(CMB)[| and the larger scale structure Q, suggest that the present universe is accelerating expansion, which needs 
a kind of mysterious dark energy with negative equation-of-state (EOS). The simplest model is by introducing the 
cosmological constant term A, which has a constant effective EOS w = —1, and drive the acceleration of the universe. 
If assuming the effective energy of the A term occupies ^ 73% of the total energy, together with 23% dark matter, 
^ 4% baryon matter and ~ 10~^ radiation component, constitute the so-called ACDM model. 

Although this simple model satisfies nearly all the cosmological observations, it still remains a phenomenological 
model. Since the major components, A and dark matter, still keep unclear for us. For the A as the candidate of dark 
energy, also suffers from the following dilemmas (see for instant ^]): First, the effective energy density is quite tiny, 
PA ~ 5.8h^ X lO^^^eV*. If we consider it as the vacuum energy of the particle physics, this energy density is 120 order 
smaller than the Planck energy scale! This is the so-called 'fine-tunning' problem. Second, the density scale of A 
keeps constant in all the stage of the universe. The observations show that the present value of the matter component 
(including dark matter and baryon components) is some one third of pA, but it varies with time as pm oc a~^. So, 
for example, at an earlier time of radiation-matter equality with redshift z 3454, pA should be a very fine tuned 
value ~ 6.3 x lO^^^pm- Otherwise, a slightly variant initial value of pA would lead to a value of the ratio Pa/ Pm 
drastically different from the observed one. This is the so-called 'coincidence' problem. In addition, the ACDM model 
also faces some observational challenges: there are mild evidences show that, the EOS of dark energy might evolve 
from w > — I in the early stage to w < — 1 in the current stage which is expected to be confirmed by the future 
sensitive observations [g*]. 

So, it is necessary to look for other candidates as the dark energy, especially the dynamical models. One possibility 
is proposing a canonical scalar field with the langrangian = '/'^/2 — V^((/>), which is dubbed as the 'quintessence' 
models (see the review f^). Similar to the inflationary field, when the potential term V{(j)) is dominant, the EOS of 
the quintessence field approaches to —1, i.e. A-like, and accounts for the observations. The most interesting is that, 
in [g], the authors introduced a kind of potential forms, such as V{(j)) = M'^~^°'(j)~°^ or V{(j)) = M*[exp(Mpi/</)) — 1], 
which have the tracker solutions, i.e. the field tracked the evolution of the background components in the early 
universe. So they address the coincidence problem, i.e. removing the need to tune initial conditions in order for the 
matter and dark energy densities to nearly coincide today. Although, this kind of models have been excluded by the 
cosmological observations, it provides the excellent possibility to naturally avoided the coincidence problems. 

However, the quintessence models also suffer from some dilemmas. The EOS of the quintessence models are in the 
range — 1 < w < 1. In order to obtain a EOS with w < —1, one always has to introduce the so-called 'phantom' 
field, which includes the non-canonical negative kinetic terms However, the 'phantom' field lacks the strong 

physical motivation, and also leads to the problem of quantum instabilities In the authors also pointed 

that, if w < — 1 keeps, the universe shall face the 'big-rip' problem. 

In addition to proposing the dynamical dark energy models, some efforts have been paid to modify the general 
relativity (GR) [l3|, which can also speed-up the universe in the recent stage. However, it should be pointed that, any 
revised GR should prepare to go through the strict solar system test , as well as to explain the various cosmological 
observations, such as the CMB temperature and polarization anisotropics power spectra. In addition, to our view, the 
current observations have not provided the strong reasons to answer: Why we should modify GR and how to modify 
it. 

In this paper, we shall propose the Yang-Mills condensate (YMC) dark energy models, where, instead of the scalar 
field, the quantum Yang-Mills field is considered as the candidate of dark energy component. Recently, the similar 
models are also discussed by a number of authors Different from the scalar field models, the Yang-Mills 

fields are the indispensable cornerstone to particle physics, gauge bosons have been observed. There is no room for 
adjusting the form of effective Yang-Mills Lagrangian as it is predicted by quantum corrections according to field 
theory. In this review, we shall firstly introduce physical motivations of the YMC dark energy models in Section II. 
In Section III, we simply introduce and discuss the Lagrangian of effective quantum Yang-Mills field. 

As the main part of this paper, in Section IV, we apply the YMC into the cosmology as the candidate of dark 
energy, and investigate the cosmic evolution of the various components, especially the evolution of dark energy. We 
find the excellent characters of this kind of models. Different from the quintessence models, both the EOS Wy > —1 
and Wy < —1 can be naturally realized. In the free YMC dark energy models, Wy 1/3 in the early stage and 
tracked the evolution of radiation component. Only in the recent stage, Wy rolls to the critical state with Wy — — 1, 
i.e. A-like, and accounts for the observations. This feature is independent of the choice of the initial condition, so 
the coincidence problem is naturally avoided. Wc also find that, if the possible interaction between YMC and dust 
components is considered, the late time attractor solution can exist, where the EOS of YMC naturally runs from 
Wj, — > 1/3 to the the phantom-like attractor state, i.e. Wy < —1. However, the total effective EOS is Wtot = —1 in the 
attractor solution, and the so-called 'big-rip' problem is also avoided. We should pointed that, these new features are 
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all independent of the interaction forms. In Section V, we calculate the statefinder and Omdiagnostics for the YMC 
dark energy models, which are helpful to differentiate the YMC models from the other dark energy models from the 
observations. 

Section VI is contributed as a summary of this paper. 

Throughout this paper, we will work with unit in which c — h — ks = 1- 



II. PHYSICAL MOTIVATION 

The introduction of the quantum effective YMC into cosmology [l7l| has been motivated by the fact that the 
SU{i) YMC has given a phenomenological description of the vacuum within hadrons confining quarks, and yet at 
the same time all the important properties of a proper quantum field are kept, such as the Lorentz invariance, the 
gauge symmetry, and the correct trace anomaly [l8| [l9|. Quarks inside a hadron would experience the existence 
of the Bag constant, B, which is equivalent to an energy density B and a pressure —B. So quarks would feel an 
energy-momentum tensor of the vacuum as T^^ — i?diag(l, — 1, — 1, — 1). This non-trivial vacuum has been formed 
mainly by the contributions from the quantum effective YMC, and from the possible interactions with quarks. Our 
thinking has been that if the vacuum inside a hadron is filled with the quantum effective YMC, what if the vacuum 
of the universe as a whole is also filled with some kind of YMC. 

Gauge fields play a very important role in, and are the indispensable cornerstone to, particle physics. All known 
fundamental interactions between particles are mediated through gauge bosons. Generally speaking, as a gauge field, 
the YMC under consideration may have interactions with other species of particles in the universe. However, unlike 
those well known interactions in QED, QCD, and the electro- weak unification, here at the moment we do not yet 
have a model for the details of the microscopic interactions between the YMC and other particle. Therefore, in this 
paper on the dark energy model based on the quantum effective YMC, we will adopt a simple description of the 
possible interactions between the YMC and other cosmic particles, in addition to the simplest model with free YMC 
component. We will investigate in these models the cosmic evolution of the universe from the radiation-dominated 
era up to the present. 



III. YANG-MILLS FIELD MODEL 

In the effective YMC dark energy model, the effective Yang-Mills field Lagrangian is given by [isj [20l. [2]|: 

C,s^\bF{\n\F/i^\-l), (1) 

where k is the renormalization scale of dimension of squared mass, F = —{l/2)F^j^F°-^^'^ plays the role of the order 
parameter of the YMC. The Callan-Symanzik coeflicient b = (UiV- 27V/)/247r2 for SU{N) with Nf being the number 
of quark fiavors. For the gauge group SU{2) considered here, one has b = ll/127r^ when the fermion's contribution 
is neglected. For the case of SU{3) the effective Lagrangian in ^ leads to a phenomenological description of the 
asymptotic freedom for quarks inside hadrons [lB|- It should be noted that the SU{2) Yang-Mills field is introduced 
here as a model for the cosmic dark energy, in may not be directly identified as the QCD fields, nor the weak- 
electromagnetic unification gauge fields. 

An explanation can be given for the form in ([1]) as an effective Lagrangian up to 1-loop quantum correct [H, [l^, . 
A classical SU{N) Yang-Mills field Lagrangian is £ = F/2gQ, where go is the bare coupling constant. As is known, 
when 1-loop quantum corrections are included, the bare coupling go will be replaced by a running one g as the following 
9o ~^ = 2/(&ln(fc^/fcn)), where k is the momentum transfer and ko is the energy scale. To build up an effective 
theory [Tsl. Il9l. [23|. one may just replace the momentum transfer k^ by the field strength F in the following manner: 
ln(fc^/A;Q) 2\n\F/eK^\, yielding equation ([1]). We would like to point out that the renormalization scale k is the 
only parameter of this effective Yang-Mills model, and its value should be determined by comparing the observations. 
In contrast to the scalar field dark energy models, the YMC Lagrangian is completely fixed by quantum corrections 
up to 1-loops, and there is no room for adjusting its functional form. This is an attractive feature of the effective YMC 
dark energy model. We should mention that, the YMC dark energy models including 2-loop and 3-loop quantum 
corrects are also discussed in the recent papers 23l[25|. It was found that, these more complicated models have the 



exactly same characters with the 1-loop models. So in this paper, we mainly focus on the 1-loop models with the 
Lagrangian in 
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IV. YMC AS DARK ENERGY 



Let us work in the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) universe, which is described by 

where t time and t time are related by cdt = adr. Considering the simplest case with only the YMC in the universe, 
which minimally coupled to the gravity, the effective action is. 



5* = 



n 



16nG 



^cff 



d^X. 



(2) 



Here, g is the determinant of the metric 5^,^. TZ is the scalar Ricci curvature, and Ccs is the effective Lagrangian 
of YMC, described by Eq. ([1]). By variation of S with respect to the metric g^'^, one obtains the Einstein equation 
Gp,i, — 87rGT^,y, where the energy- momentum tensor of YMC is given by 



(3) 



The dielectric constant is defined by e = 2d£cff/dF. In the one-loop order, it is given by 

e = bln\F/K^\ . 



(4) 



This energy-momentum tensor (|3]) is the sum of three different energy-momentum tensors of the vectors, T^^ = 
J2a '^"'^T^y, neither of which is of prefect-fluid form. Here we assume the gauge fields are only the functions of time 
t, and — ^(^^^"(t) (here CTq are the Pauli's matrices) are given by = and A° = 5fA{t). Thus, we will find 
that, the total energy-momentum tensor T^^, is homogeneous and isotropic. 
Define the Yang-Mills field tensor as usual: 



(5) 



where /"^'^ is the structure constant of gauge group and /"^'^ = e"''^ for the SU{2) case. This tensor can be written 
in the form with the electric and magnetic field as 



/ El E2 E3 

—El B3 —B2 

-E2 — -B3 Bi 

\-E3 B2 -Si 



(6) 



It can be easily found that Ef = E2 = E^, and Bf — B2 — B^. Thus F has a simple form with F = E^ — B^, where 
E'^ — J2^=i ^^'^ = X]^=i ^i- this case, each component of the energy-momentum tensor is 



Wr;^ = lie-b)iB'-E')6l + lEX, 

^T} ^ \{e- b)iB^ - E^)d] + ^E^S}S] - ^B^S}{1 - S^). 

D 00 



(7) 
(8) 



Although this tensor is not isotropic, its value along the j = a direction is different from the ones along the directions 
perpendicular to it. However, the total energy-momentum tensor = X]a=i ^""^T^j^u has isotropic stresses, and the 
corresponding energy density and pressure are given by 



py = \e{E^ + B^) + \b{E^ B% py = le{E' + B^) \b{E^ B^). 



(9) 



In this paper, for simplicity, we only discuss the pure 'electric' case, F = E^ . This is a typical consideration, since in 
the expanding universe, a given 'magnetic' component of Yang-Mills field decreases quite rapidly, and the Yang-Mills 
field becomes the 'electric' type [2^. The energy density and pressure of YMC are reduced to 



E^ 



Py 



{e + b), Py 



E^ 
2 
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It is convenient to introduce a dimensionless quantity y = e/h = ln|_F/K^|. So the quantities py and py can be 
rewritten as 

Py = \hn\y + l)ey, ^ lhn\y - i)ey . (10) 

One sees that, to ensure that energy density py be positive in any physically viable model, the allowance for the 
quantity y should he y > -1, i.e. F > n^/e ~ 0.368^2. The EOS of the YMC is 

^ y-^ 

' 3y + 3 ^ ^ 

Before setting up a cosmological model, the EOS Wy itself as a function of F is interesting. At the state with F = k'^, 
which is called critical point, one has y = and the YMC has an EOS of the cosmological constant with Wy ~ — 1. 
Around this critical point, F < gives y < and w < —1, and F > gives y > and w > —1. So in the YMC 
dark energy models, EOS with w > —1 and w < — 1 all can be naturally realized. On the other hand, in the high 
energy scale with F k^, one finds that the YMC exhibits an EOS of radiation with Wy = 1/3. In the follows, we 
will detailed show that the YMC was evolving from the state with Wy = 1/3 to = —1 (or even to Wy < —1) in 
the expanding universe. This is the main context in this paper. In addition, the characteristic statefinder parameters 
and the perturbations of the YMC dark energy is also discussed in [13, HI] [1^ , which are helpful to distinguish the 
YMC model from other dark energy models. 

As is known, an effective theory is a simple representation for an interacting quantum system of many degrees of 
freedom at and around its respective low energies. Commonly, it applies only in low energies. However, it is interesting 
to note that the YMC model as an effective theory intrinsically incorporates the appropriate states for both high and 
low temperature. As has been shown, the same expression in (|lip simultaneously gives Py ^ ~py at low energies, 
and Py Py/3 at high energies. Therefore, the model of effective YMC can be used even at higher energies than the 
renormalization scale k. 



A. Free YMC models 



Let us discuss the cosmological model, which filled with three kinds of major energy components, the dark energy, the 
matter, including both baryons and dark matter, and the radiation. Here, the dark energy component is represented 
by YMC, and the matter component is simply described by a non-relativistic dust with negligible pressure, and the 
radiation component consists of the photons and possibly other particles, such as the neutrino, if they are massless. 

Since the universe is assumed to be flat, the sum of the fraction densities is fly + f2m + — 1, where the fractional 
energy densities are Qy = Py/ptot, = Pm/ptot, = Pr/ptot, and the total energy density is ptot = Py + Pm + Pr- 
The overall expansion of the universe is determined by the Friedmann equations: 



a J 3 
& AttG 



8ttG , , , , 

[Py+Pm+Prj, (12) 



a 3 



{Py + 3py + Pm + Pr +SPr), (13) 



where and what follows, the superscript dot denote the d/dt. These three components of energy contribute to the 
source on the right-hand side of the equations. We should notice that, Pr — Pr/3- In this subsection, we shall 
assume there is no interaction between these three components. The dynamical evolutions of the three components 
are determined by their equations of motion, which can be written as equations of the conservation of energy. 



Py+3%y+Py)^0, (l4) 

a 

Pm + S-pm = 0, (15) 

a 

Pr+3-{pr+Pr)^0. (16) 

a 



As is known, Eq. is not independent and can be derived from Eqs. (HH), (fTS]) . (|16p and p^ . From Eqs. (fT5 
and (jl6p . we can obtain the evolutions of the matter and radiation components, p„i cx a^^ and p,n oc a^^. 

Now, let us focus on the evolution of YMC. Inserting pO|) into (fT4|) . we can obtain the following simple relation. 
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FIG. 1: In the free YMC dark energy models, the evolution of the YMC energy density (left panels) and EOS (right panels) 
for the models with different initial conditions. 



where the coefficient of proportionality in the above depends on the initial condition. At the early stages, a ^ 0, 
Eq. (fTT]) leads to y » 1, and Eq. pT|) gives Wy 1/3, so the YMC behaves as the radiation component. With 
the expansion of the universe, the value of y runs to the critical state of y = 0, and the EOS goes to Wy — ~1. So, 
in the late stage of the universe, the YMC behaves as the cosmological constant. This is one of the most attractive 
character of the YMC models. Around the critical point with y — 0, Eq. (fT7|) yields y cx a~^, and the EOS of YMC is 
+ 1 ~ 4y/3 cx a~^. The YMC can achieve the states of Wy > —1 and Wy < —1, but it can not cross over —1, just 
like in the scalar models [30| . 

We should notice that, the relation in (fT7|) can also be derived by the effective Yang-Mills equations. By variation 
of the action S with respect to gauge field A^, we can obtain the Yang-Mills equations, which are 

d^,{a^e F"""") + f"'^'' A^ia'^e F"'^") = 0. (18) 

The V = component of which is an identity, and the v — 1,2,3 spatial components can be reduced to dria^^E) = 0. 
At the critical point (e = 0), this equation is an identity. When e 7^ 0, this equation has an exact solution as pT)) [2l| . 
Now, let us fix the value of k, the only parameter in the model. At the present time, the YMC energy density is 

py = — (y + 1)-— , (19) 

and, as the dark energy, it should be ^lyptot, where the present total energy density in the universe ptot ~ 8.099 x 
10~^^eV''. We choose fly — 0.73 as has been observed, yielding 

K = 3.57/1 X lO^^eV^ (20) 

This energy scale is low compared to typical energy scales in particle physics, such as the QCD and the weak- 
electromagnetic unification. This is the reason, why the SU{2) Yang-MilUs field introduced in this paper cannot be 
directly identified as the the QCD gluon fields, nor the weak-electromagnetic unification gauge fields. 

To be more specific about how the YMC evolves in the expanding universe, we look at an early stage when the Big 
Bang nucleosynthesis (BBN) processes occur around a redshift z ~ 10^° with an energy scale ~ IMeV. To see how 
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FIG. 2: Left panel: The present EOS of the YMC wo depends on the initial value rj. Right panel: The transition redshift zt 
depends on the initial value rb- In both panels, the shadow region has been excluded by the observation of BBN. 

the evolution of py depends the the initial condition, we introduce the ratio of energies of the two components 



where pr is the radiation energy density. We consider rt, < 1, i.e. the YMC is subdominant to the radiation component 
initially. Of course, the YMC evolves differently for different initial values of r^. Nevertheless, we will see that, as the 
result of evolution, the present universe is always dominated by the YMC ^ 0.73 for a very wide range of initial 
values ffc. 

Now we use the exact solution p7)) to plot the evolution of py as a function of the redshift z in Fig. [l](left panel). 
As specific examples, here we take n = 10-2, and n = lO'^o. in comparison, also plotted are the energy densities 
of radiation, and of matter. It is seen that, in the early stages, py decreases as Py cx a^"'. So the YMC density is 
subdominant and tracks the radiation, a scaling solution. The corresponding EOS of Yang-Mills field is Wy ~ 1/3 
shown in Fig. [1] (right panel). At late stages, with the expansion of the universe, a oo, y decreases to nearly zero, 
and lUj, — !■ — 1 asymptotically. Moreover, this asymptotic region is arrived at some redshift z before the present time, 
and this z has different values for different initial values of rf,. For smaller rf,, the transition redshift is larger (seen in 
Fig. [T]), and the transition happens earlier. Once the asymptotic region is achieved, the density of the YMC levels off 
and remain a constant forever, like a cosmological constant. We have also checked that the present value fly ^ 0.73 is 
also the outcome of the cosmic evolution for any value of rf, in the very wide range (10^^*^, 10^^). So the coincidence 
problem do not exist in the YMC model. 

The present EOS of the YMC wa is nearly —1. Fig. [2] (left panel) plots the dependence of the present EOS wq on 
the initial condition rt- The function log]^Q(rf,) versus log2Q(u;o + 1) is nearly linear: a smaller rt, leads to a smaller 
Wq. For a value ri, — 10"^, one has wq — —0.99. For a value rf, = 10"^°, wq would be —1 accurately up to one in 
10^^. Therefore, at present the YMC is very similar to the cosmological constant. 

The solution in Eg.pT)) can converted into the following form 




(21) 



Pr z=10io 




(22) 



where yo is the value of y at z = 0, depending on the initial value rf,. For a fixed yo, this formula tells a one-one 
relation between the EOS (through y) and the corresponding redshift z. As is seen from Fig. [TJ the transition of lo 



from 1/3 to —1 occurs during a finite period of time, instead of instantly. To characterize the time of transition, we 
use Zt to denote the redshift when w = —1/3, i.e. y = 1, as given by Eg. (fTTj) . This is, in fact, the time when the 
strong energy condition begins to be violated, i.e., py + ipy < 0. Then 




1 - 2/0 



- 1. (23) 



Therefore, this gives a function zt — zt (?"(,). Fig. [2] (right panel) shows how the transition redshift zt depends on 
the ratio rf,. Interestingly, this transition can occur before, or after the radiation-matter equality {zeq = 3454). This 
feature is different from the tracked quintessence models in which transition occurs during the matter dominated era 
Q. A larger rt leads to a smaller Zf. For example, rt — 10^^ leads to zt ~ 12.4 ^ z^q, and the transition occurs in 
the matter dominated stage, and 10~^° leads to zt ~ 5.0 x 10^ 3> z^q, and the transition occurs in the radiation 
dominated stage. 

The value of cannot be chosen to arbitrarily large. In fact, there is a constraint from the observation result 
of the BBN. As is known, the presence of dark energy during nucleosynthesis epoch will speedup the expansion, 
enhancing the effective species of neutrinos [3l| . The latest analysis gives a constrain on the extra neutrino species 
SN^ = — 3 < 1.60 [3l|. Here in our model, the dark energy is played by the YM field. By a similar analysis, the 
ratio Tfc is related to SN^ through rf, — ^^^^j^. This leads to an upper limit < 0.26, the present EoS < —0.94 
by Fig. [2] (left panel), and the transition redshift zt > 5.8 by Fig. [2] (right panel). The range of initial conditions 
Tf, e (10~^", 10^^) that we have taken satisfies this constraint. 

B. Coupled YMC models 

In this subsection, we shall generalize the original YMC dark energy model to include the interaction between the 
YMC and dust matter. We should mention that, the possible interaction between YMC and radiation component may 
also exist, which has been discussed in [33|- (The similar models on the scalar field dark energy have been discussed 
by a number of authors (see [s^ for instant)). In this section, we assume the YMC dark energy and background 
matter interact through an interaction term Q. Thus the equations of the conservation of energy in (|14p and (jl5p 
should be changed into 

f>y + 3\y+Py)^-Q, (24) 

Pm+3^Pm^Q, (25) 

and the equation for radiation in pB]) is still held. The sum of Eqs. and (pS)) guarantees that the total energy 
of YMC and dust matter is still conserved. It is worth noting that the free Yang-Mills equation in ([T5|) is not 
satisfied when Q ^ 0. In the natural unit, the interaction term Q has the dimension of [energy]^. The coupling Q is 
phenomenological, and their specific forms will be addressed later. When Q > 0, the YMC transfers energy into the 
matter, and this could be implemented, for instance, by the processes with the YMC decaying into pairs of matter 
particles. On the other hand, when Q < 0, the matter transfers energy into the YMC. Note that, once Q is introduced 
as above, it will bring another new parameters in the models, in addition the free parameter k. 

We introduce the following dimensionless variables, rescaled by the critical energy density 6k^/2 of the YMC, 

_ 2pm , 2(5 

where / is the function of x and y. By the help of the definition of y, the evolution equations and (P5|) can be 
rewritten as a dynamical system, i.e. 

/ _ 4y _ f{x,y) 
y - 2 + y (2 + y)ev' ^''^ 
x' = -3x + /(x,2/). (28) 

Here, a prime denotes derivative with respect to the so-called e-folding time N = In a. The fractional energy densities 
of dark energy and background matter are given by 



{l + y)ey + x' {l+y)ev+x' 



= , ^.^:„ , ^ , and n„, = ^ ^ ^ . (29) 
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Before discussing the specific form of the interaction term Q, let us first investigate the general feature of this 
dynamical system, described by ([27|) and ([28| . We can obtain the critical point {yc,Xc) of the autonomous system by 
imposing the conditions y'^ = x'^ = 0. From the equations ([27|) and (pSj) . we obtain that the critical state satisfies the 
following simple relations 

Sxc = f{xc,yc), (30) 
3xc = -42/ee^=. (31) 

So we can obtain the critical state (j/c, Xc) by solving these two equations. In order to study the stability of the critical 
point, we substitute linear perturbations y ^ yc + Sy and x ^ Xc + Sx about the critical point into dynamical system 
equations (|27p and (|28p and linearize them. Thus, two independent evolutive equations are derived. 



dx' J V / V J^y /,x - 3 y \^ da; 

where R^y = dR/dy at {y,x) — (yc,Xc)- The definitions of R ^, /,y, /.x and G^y are similar. The functions G and i? 
are defined by 

G^Giy)=-^, R.Rix,y)- ^^^''^ 



y (2 + y)e^' 

which are used for the simplification of the notation. The two eigenvalues of the coefficient matrix Af determine 
the stability of the corresponding critical point. The critical point is an attractor solution, which is stable, only if 
both the these two eigenvalues are negative (stable node), or real parts of these two eigenvalues are negative and the 
determinant of the matrix M is negative (stable spiral). 

Here we discuss some general features of the attractor solutions, regardless the special form of the interaction term 
Q. From the expression ((3T|l . we find that Xc = — (4/3)?/ce^''. Substitute this into the formula ([29|l . one obtains 



{yc + l)e^- ^ 3 + 3ye 
{yc + l)ey'= + Xc i-yc 

Since < fij, < 1, this formula follows a constraint of the critical point 

- 1 < J/c < 0. (33) 
From the formulae (fTT|) and (|32p . we find a simple, but interesting relation, 

nyWy = -1. (34) 

This relation is held for all attractor solutions, independent of the special form of the interaction. Since the value of 
fij, is smaller than or equal to unity in the attractor solution, we obtain that 



Uy 



< -1- (35) 



This means that, the EOS of the YMC dark energy must be smaller than or equal to —1, i.e. phantom-like or A-like. 
Since in the early universe, the value of the order parameter of the YM field F is much larger than that of k^, i.e. 
J/ 3> 1, the YM field is a kind of radiation component. However, in the late attractor solution, the dark energy 
is phantom-like or A-like. So the phantom divide must be crossed in the former case, which is different from the 
interacting quintessence models. 

It is interesting to investigate the total EOS of the YMC and dust matter, which determines the finial fate of the 
universe, when the radiation becomes negligible. The total EOS is defined by 

Plot Py Pm ^ for^ 

Wtot = = ; = >'lyWy, (36) 

Ptot Py + Prn 



where pm = is used. From the relation p4[) . we obtain that, in the attractor solution, 

Wtot = -1- (37) 

This result is also independent of the specific form of the interaction. So the universe is an exact de Sitter expansion, 
and the cosmic big rip is naturally avoided, although the YMC dark energy can be phantom-like. 



Now, let us discuss the evolution of the various components in the universe for some specific interaction models. 
In this paper, we shall focus on the following three phenomenological interaction models: Q oc Hpy, Q oc Hp^j, and 
Q oc H{py -\- Pm): separately. Some other phenomenological models have also been discussed in the papers 33, 34 1. 
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1. Q QC Hpy 



In this case, we can write the function f(x,y) as the following form, f{x,y) — a{l + y)e^ ■ Of course, when a = Q. 
the system returns to the model with free YMC dark energy. Here, we consider the simplest case with a being a 
non-zero dimensionless constant. Thus, the dynamical equations in (P7)) and pS)) becomes 

Ay a[l + y) 

y - ^7^. 77;-, — T' ^"^^i 

2 + y {2 + y) 

x' = -3x + a{l + y)ey . (39) 

Obviously, the evolution of dust matter and YMC are influenced by the interaction by the function f{x,y). We 
expect, when the fraction density of YMC was sub-dominant in the universe, the effect on the dust is small. Only in 
the latest stage of the universe, where the YMC dark energy dominates the evolution of the universe, the effect of 
interaction on the dust becomes important. 

The critical point (j/c, Xc) is obtained by imposing the condition y'^ — x'^ — 0, which are 

yc = -j^, xc = -^ey', (40) 

4: + a 3 

and the fractional energy density and the EOS of the YMC at this critical point are 

^y = -— = (41) 
Wy a + 6 

The constraint < fij, < 1 requires that a > (note, we have set a ^ throughout the discussion). In this condition, 
we find that Wj, < — 1 is satisfied. In order to keep this critical point being stable, i.e. (jjc, Xc) is the attractor solution, 
a constraint on a can be derived: a > —8, which is auto-satisfied when a > is required. So the attractor solution 
requires that the constraint on the coefficient being 

a > 0, (42) 

which follows that the EOS of YMC in this attractor solution must be negative. If we still require that, the fraction 
density of YMC is larger than the value at the present stage, i.e. fly > 0.73, the constraint on a becomes much tighter 
< a < 1.11. 

In order to have a much clear picture for this system, let us study the evolution of the various components in the 
universe by adopting a = 0.5. We still choose the initial condition at the BBN stage, i.e. z — 10^°, where the value 
of X is chosen as Xi = 1.8 x 10~^^ to keep the present value of fly being fij, — 0.73. For the YMC, we consider the 
following two choices as the initial condition: i.e. yi = 60 and = 20. In the former case, we have VLyi/fl^i — 0.16, 
and in the latter case, we have ^yi/^mi = 3 x 10^^^. Although, the difference between these two flyi/Q.mi cases is 
larger than 17 orders, we will show that, these two models follow the similar present state of universe. 

We should also fix the value of k. In order to keep the present total energy density in the universe being ptot — 
8.099/i2 X lO-^eV*, we find the n = 4.00/i x lO'^eV^ in both initial condiitons. If the Hubble constant h = 0.72 is 
adopted, we get k = 2.88 x 10~^eV^. Again, we find the value of k is much smaller than the typical energy scale in 
the particle physics. 

In Fig. [3] (left panel), we plot the evolution of various components in the universe. Note that, the evolution line for 
dust component covered each other in both initial cases, as well as the radiation component. Similar to the free YMC 
models (shown in the left panel in Fig. [T|), in the early stage of the universe, the YMC was tracking the radiation, 
and it transfered to the attractor solution in the later stage. The smaller yi induces the earlier transfer. In the right 
panel of Fig. [31 we plot the evolution of EOS of YMC Wy and the total EOS Wtot- Note that, the evolution line Wtot 
covered each other in both initial cases. As expected, in the both initial conditions, in the early stage, Wy 1/3. 
In the late stage, Wy runs to the attractor solution with Wy = —{a + 3)/3 with a — 0.5, i.e. Wy —^ —7/6. In the 
intermedial stage, the EOS Wy must cross the line with Wy = —1. For the total EOS wtot, in the early stage, Wtot 0, 
where dust component is dominant than YMC. However, in the late stage, Wtot ^ ^1, which is consistent with the 
previous discussions. 
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log„(1+z) log„(1+z) 



FIG. 3: In the coupled YMC dark energy models, the evolution of the YMC energy density (left panels) and EOS (right panels) 
for the models with different initial conditions. 



•2,. Q QC Hpm 

Now, let us turn to the interaction cases with Q cx Hpm, which is equivalent to the form f{x,y) — (3x. The 
dynamical equations in ([27ll and ([28|) becomes 

, Ay l3x , 

x' = {P-3)x. (44) 

If we consider the simplest, where is a constant, the equation (|44|) follows that x oc a'^^^. From the definition of x, 
we derive that pm oc a'^"^. When (3 — 0, the model returns to the free YMC cases, and pm oc as usual. However, 
when /3 7^ 0, the evolution of the dust component is changed, which is conflicted with the evolution of dust in the 
standard hot big-bang model. So it is dangerous to consider this kind of interaction term in the early stage of the 
universe. 

However, it is allowed to consider the form of Q oc Hp,n as a kind of phenomenological model in the late stage of 
the universe. The critical point (yc,Xc) of (j43p and (j44p is obtained by imposing the condition y'c = x'^ — 0, we find 
these is no solution at all. So we conclude that, it is impossible to obtain an attractor solution for this kind of system. 



3. Q cx H{py + p™.) 

In the end, let us discuss the another kind of phenomenological model, where Q oc H{py + p^) is satisfied. This is 
equivalent to set the form /(x, y) — 'y[(y + l)e^ + x]. We consider the simplest case, where 7 is a non-zero dimensionless 
constant. In the dark energy dominant stage, this system returns to the first case with Q oc Hpy, and in the dust 
dominant stage, it returns to the second case with Q oc Hp„i. Similar to the second case, if we directly apply this 
system to the early universe, the evolution of dust would be changed, which is conflicted with the prediction of the 
standard hot big-bang models. 

In this paper, we also consider this interaction form as a kind of phenomenological model in the late stage of the 
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FIG. 4: Left panel: The r — s diagram of the YMC models; Middle panel: The r — g diagram of the YMC models; The right 
panel: Om{z) diagnostic in the YMC models. 



universe. The dynamical equations in ([27|) and ((28|) becomes 

y 2 + v {2 + y)ey ' ^^"^ 

x' = -3a; + 7[(?/ + l)e^ + a;]. (46) 

From the equations pO|) and (|3T|) . we obtain the critical point 

2/c = ^, -c=-%e^^ (47) 
The fractional energy density and the EOS of the YMC at this critical point are 

The constraint of < fij, < 1 requires that < 7 < 3. In order to keep this critical point being stable, i.e. (?/c, Xc) is 
the attractor solution, another constraint of 7 can be derived: 7 > —120/31, which is auto-satisfied when < 7 < 3 
is required. So the attractor solution requires that the constraint on the coefficient being 

< 7 < 3, (49) 

which follows that, the EOS of YMC in this attractor solution must be negative. If we still require that, the fraction 
density of YMC is larger than the value at the present stage, i.e. Vly > 0.73, the constraint on a becomes much tighter 
< 7 < 0.81. 



V. STATEFINDER AND Om DIAGNOSIS IN THE YMC MODELS 



In this section, we shall present the way to discriminate between the YMC models and the other dark energy 
models. In the previous works [3^[3^, the authors suggested to use the so-called "statefinder" {r,s} pair and Om 
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diagnostics. Here, we shall also apply the statefinder diagnosis into the YMC models. The similar discussion can be 
found in the previous works [l^l [lH [ll] [s^ • 

The statefinder diagnostic pair {r, s} is defined as 

(50) 



~ai73' ~3(g-l/2)' 

where q = —a/iaH^) is the deceleration parameters, r forms the next step in the hierarchy of geometrical cosmological 
parameters beyond H and q, and s is a linear combination of r and q. Apparently, the statefinder parameters depends 
only on scale factor a and its derivatives, and thus it is a geometrical diagnostic. These parameters can also been 
expressed in terms of ptot and ptot as follows 

1 /, , 'iPtot\ , 9{ptot + Ptot)ptot Ptot + Ptot Ptot , 

9=0 1 + , r = l + — , s = — . (51) 

^ \ Ptot ) IptotPtot Ptot Ptot 

By using Eqs. and (|25p . we plot the trajectories of (r, s) and of (r, q) in Fig. |3]from the redshift z = 10, where 
we have considered two cases with Q = and Q = O.bHpy. The arrows along the curves indicate the direction of 
evolution. In both cases, we have chosen the initial condition at the redshift z — 3454, where the densities of matter 
and radiation equate to each other. We choose the initial condition to make the fij, = 0.01 at this high redshift. 

From this figure, we find that, the model with different interaction terms have the difference in both (r, s) and (r, q) 
trajectories. However, in the low redshift (z < 10), especially when z 0, the trajectories in both models become 
quite close to those of the ACDM model. However, we should also mention that, the trajectories of the statefinder 
shows in this figure are quite different from the other dark energy models, such as the decaying vacuum model, the 
Quintessence models, the K-essence models (see for instant [33| ). In order to break the degeneracy between the YMC 
models and the ACDM models, let us consider the Om diagnostic, which is defined as 

Om{x) = '^yi^^ (52) 

where x = (1 + z), and h{x) = H{x)/Hq. Thus Om involves only the first derivative of the scale factor through the 
Hubble parameter and is easier to reconstruct from the observational data. For the ACDM model, it is simple, i.e. 
Om{x) = r^m, independent of the redshift. For YMC model, we plot Om as a function of z in the FigH] (right panel). 
It is interesting to find that, for either model, Om is nearly a constant at the low redshift. For the model without 
interaction, as expected, Om ~ 0.27, very close to that of ACDM models. However, for the model with Q ~ 0.5Hpy, 
Om = 0.15, which is helpful to differentiate the coupled YMC models from the ACDM model. 



VI. CONCLUSION 



In order to answer the observed accelerating expansion of the universe in the present stage, we introduce the quantum 
Yang-Mills condensate dark energy models. Different from the general scalar field models, the quantum Yang-Mills 
fields are the indispensable cornerstone to particles, and the Lagrangian of the Yang-Mills field is predicted by quantum 
corrections according to field theory. 

In this paper, we review the main characters of the YMC dark energy models, where both free YMC model and 
possible coupled YMC models are considered. In all these models, the EOS of YMC is close to Wy = 1/3 in the early 
stage, and tracked the evolution of radiation components. In the late stage, the value of Wy runs to the attractor 
solution with Wy = —1 for the free YMC model, or with Wy < ~1 for the coupled YMC models. This naturally 
explains the observations. The present state of the universe is independent of the initial state of Yang-Mills field, 
and the coincidence problem is naturally avoided. In the coupled YMC models, not only the state of Wy < —1 can 
be naturally realized, mildly suggested by observations, but also the 'big-rip' problem is auto-avoided. The most 
importance we should mention is that, as a dynamic dark energy model motivated by quantum effective YM field 
Lagrangian, all the YM models based upon 1-loop, 2-loop, and 3-loop quantum corrections, respectively, have similar 
dynamic behaviors. That is, the main properties of YM model for dark energy remain stable when the number of 
loops for quantum corrections increases up to 3-loop. 

However, in all the models, the fine-tunning problem still exists, i.e. the value of k ~ 10~^eV^, which is too low 
comparing with the typical scales in particle physics. We should notice that, this problem exists in nearly all the dark 
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energy models, which may suggest a new physics at this low energy scale. 
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